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Abstract. We consider an extended Korteweg-de Vries (eKdV) equation, the usual 
Korteweg-de Vries equation with inclusion of an additional cubic nonlinearity. We 
investigate the statistical behaviour of flat-top solitary waves described by an eKdV 
equation in the presence of weak dissipative disorder in the linear growth/damping 
term. With the weak disorder in the system, the amplitude of solitary wave randomly 
fluctuates during evolution. We demonstrate numerically that the probability density 
function of a solitary wave parameter k which characterizes the soliton amplitude 
exhibits loglognormal divergence near the maximum possible k value. 
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1. Introduction 

It is well known that the Korteweg-de Vries (KdV) equation governs the propagation of 
shallow water waves of moderately small amplitude, where a balance between quadratic 
nonlinearity and linear dispersion results in soliton solutions. For the waves of large 
amplitudes or under certain circumstances in stratified fluids, however, it was found 
that an additional cubic nonlinearity becomes crucial, leading to the extended KdV 
(eKdV) equation [U EJ [3] . The eKdV equation, also known as the Gardner equation, 
appears as a governing equation for long interfacial waves in two-layer system [3J HI [5] 
as well as for oceanic stratification in shear flow [HI [71 El E] • 

Both the KdV and eKdV equations are exactly integrable. The integrability induces 
that the eKdV equation possesses conventional soliton solutions of small amplitudes 
similar to those of the KdV equation. In addition, this extended equation exhibits 
a remarkable feature that distinguishes itself from the KdV equation, namely, the 
emergence of large amplitude wide solitons, called flat-top solitons. Such wide solitons 
also appear as solutions for other related nonlinear evolution equations led by a balance 
between dispersion and nonlinearity, such as high-order Nonlinear Schrodinger equations 
(NLS) or cubic-quintic complex Ginzburg-Landau equation. Consequently, flat-top 
solitary waves draw much attention in various areas of physics including fluid mechanics 
[H QUI HH Q21 03! , nonlinear optics [H EEl EE] , and plasma physics pH EE] . 

Recently, we investigated the effects of weak dissipative disorder on flat-top solitary 
waves in cubic-quintic nonlinear Schrodinger equation (CQNLSE) and the derivative 
CQNLSE (DCQNLSE) [19], [20]. In particular, two most common types of disorder, 
disorders in the linear and cubic nonlinear gain/loss coefficients, are considered, which 
lead to the random variation of the solitary wave parameters including amplitude 
and group velocity. In this study, we showed numerically and analytically that the 
probability density function (PDF) of the soliton amplitude exhibits loglognormal 
divergence near the maximum possible amplitude value. This phenomenon stems 
from the fact that the solitary wave obtains a typical table-top shape when its 
amplitude approaches the maximum possible value. Thus, our anticipation was that the 
loglognormal divergence of amplitude can be generally associated with emergence of flat- 
top solitary waves. In the current paper, we corroborate this generality by showing that 
a perturbed eKdV equation whose structure differs from the CQNLSE and DCQNLSE, 
demonstrates similar statistical behaviour, namely, loglognormal divergence. 

The eKdV equation is not of the NLS type which was considered in our previous 
studies. On the other hand, it is ubiquitous and belongs to one of the families 
of integrable nonlinear partial differential equations. Therefore, it is important to 
examine if the loglognormal divergence of solitary wave parameters can be also found 
in a perturbed eKdV equation in the presence of weak disorder. In particular, we 
focus on the case when the dissipative disorder appears in the linear growth/damping 
term and its intensity is weak so that the solitary waves can evolve without severe 
distortion. This type of disorder can emerge quite commonly in systems that involve 
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nonlinear wave equations. Indeed, for the case of water waves, as the depth of channel 
gradually increases or decreases, the evolution of waves can be described by a perturbed 
KdV equation, where a linear term proportional to the wave envelope is incorporated 
[21~1 [22| |23"1 [23] . In the context of nonlinear optics, the random variations in the gain of 
amplifiers which are positioned to compensate the loss can lead to disorder in the linear 
gain coefficient [25J. Such disorder also appears in massive multichannel transmission 
systems due to the interplay of Raman cross talk and bit pattern randomness [26l [27] . 

Considering a perturbed eKdV equation where a random disorder appears in the 
linear growth/damping term, the amplitude of the solitary wave undergoes random 
fluctuations during evolution. Thus, we conduct Monte Carlo simulation to achieve the 
PDF of a parameter characterizing the soliton amplitude and verify its loglognormal 
divergence. This finding in turn, concludes that the loglognormal divergence of the 
amplitude PDF found in [191 EQ] is not restricted to the solitary waves of NLS type 
equations. We also note that theoretical analysis of the perturbed KdV type equations 
is an extremely challenging task due to the substantial effects of radiation. More 
specifically, the linear perturbation induces a shelf consisting of radiative modes directly 
behind the solitary wave. While the shelf has a slowly varying small amplitude, its range 
extends with time, which varies at the rate of an order one J3j [211 12H] • This phenomenon 
brings most of difficulties associated with the theoretical analysis of perturbed KdV type 
equations, and many questions still remain open despite various theoretical methods 
available [211 1231 [291 EO] ■ Although we expect that the adiabatic perturbation technique 
employed for the models in our previous study [20] can be an appropriate tool to deal 
with the underlying problem, it calls for more extensive perturbative calculations in 
order to incorporate the full impact of shelf. We therefore defer the complete theoretical 
analysis to a future publication. 

The material in this paper is organized as follows. In Sec. 2, we briefly describe the 
evolution of solitary waves of an eKdV equation in the presence of disorder in the linear 
growth/damping. In Sec. 3, we present the results of direct numerical simulations. 
Finally, in Sec. 4, we summarize our main results. 

2. Extended Korteweg-de Vries equation with disorder in the linear 
growth/damping coefficient 

We consider the evolution of solitary waves described by an eKdV equation with disorder 
in the linear growth/damping coefficient, 

d t u + 6w(l — t n u)d z u + d^u = e£(t)u. (1) 

In the context of internal waves, u represents the amplitude of the wave (or the interfacial 
displacement), z is the horizontal coordinate, and t is time. The right hand side term 
e£(t)u is responsible for the disorder effects of the linear growth/damping, and e n is 
the cubic nonlinear coefficient. We assume that the disorder £(£) is zero in average and 
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Figure 1. Solitary wave solutions u s (z,t) at t = for different n values. From 
the lowest to the highest amplitude solutions, each solution corresponds to k — 
0.5, 0.624, 0.624 999, 0.624 999 99. 



short correlated in time, i.e., 

<£(*)> = o, mm) = D6(t-t'), (2) 

where D is the disorder intensity. 

When e = 0, we obtain an unperturbed eKdV equation whose soliton solutions are 
given by 

An 2 

Us[Z ' } ~ (l-K 2 /«2 l ) 1 /2 CO sh(2 a ;) + l' [) 



where x = k(z — 4k 2 £), n m = 0.5/y^. Using the solution form we find that the 
parameter k characterizes the soliton amplitude and group velocity with the relation, 



K = y/rj/2- WOKJ, (4) 
where rj represents the soliton amplitude. These solitons are limited in amplitude 
and speed, namely, as k becomes close to its maximum possible value K m , the soliton 
forms the flat-top shape and the limiting flat-top soliton corresponds to the maximum 
amplitude and speed. Figure [1] illustrates the solitary wave solutions (EJ), for different 
k values that range from 0.5 to 0.624 999 99, where we take e n = 0.64 corresponding to 
n m = 0.625. 

In the presence of perturbation, i.e., with nonzero e, a shelf is generated in the 
lee of the solitary wave and extends its range with time. Note that this interesting 
feature imposes major difficulties on theoretical analysis of Here, we illustrate the 
emergence of shelf by numerically integrating ([1]) for a given disorder realization with 
D = 3 and e = 0.09. For the initial condition in the form of ([3]) with k = 0.5 and 
e n = 0.64, figure [2] demonstrates the solution u(z,t) at t = 10 which consists of a 
solitary wave and a shelf. 
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Figure 2. The solution u(z,t) at t = 10 for a given disorder realization £(t) with 
D = 3 and e = 0.09. The initial condition u(z,0) is in the form of ((3|) with k = 0.5 
and e„ = 0.64. 



3. Numerical simulation 

We conduct Monte carlo simulations for (Q]) with about 100 000 disorder realizations. 
Our initial condition is in the form of the solitary wave solution u s (z,0) with k = 0.5 
and e n = 0.64 which corresponds to K m = 0.625. We also assume that the disorder 
intensity is weak, in particular, D = 3 and e = 0.09. The numerical simulation is 
carried out until the solitary wave reaches a time tf — 10, where the disorder strength 
becomes De 2 tf = 0.243. We integrate (TjQ) by employing a fourth-order split-step 
method which combines the fast Fourier transform and a fourth order Runge-Kutta 
scheme for the linear and the nonlinear parts of the equation, respectively [31]. We also 
introduce artificial damping near the boundaries of computational domain in order to 
avoid numerical artifacts resulting from the radiation emission and the use of periodic 
boundary conditions [32]. The size of computational domain is taken large enough, 
— 100 < z < 100, so that the damping layer does not affect the major portion of solitary 
wave dynamics. The descritized time and spatial steps are taken as 5t = 0.001 and 
5z = 0.024, respectively. 

After retrieving the shape of solitary wave at the end of the evolution, we calculate 
the value of k based on (jlj). Repeating this procedure for independent realizations of 
disorder, we achieve the PDF of k shown in figure [31 The result of numerical simulation 
clearly demonstrates that the PDF attains a divergence near K m . 

We now explain how to verify the divergence observed in figure [3] is indeed 
loglognormal. First, it should be mentioned that a preliminary study has obtained 
an analytic form of k-PDF at the first order adiabatic perturbation (see e.g., [20] for 
the details of the adiabatic perturbation technique). Due to the emergence of shelf 
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and its nature, however, a higher order perturbation theory is necessary for more 
exact description of the PDF. Nevertheless, the result of the first order perturbation 
calculation suggests that the analytic form of k-PDF denoted by F(k) approximately 
follows (see e.g., (7) in [20J), 

^ ^ exp{ - A 2 In 2 [Aiarctanh(/t/ K m )]} 
K m (l - /t 2 //4Jarctanh(fi;/K m ) 

for < k < K m and F(k) = elsewhere. Here, Ai, A 2 , A3 are some constants related to 
the parameters D, e, and the total evolution time tf. This expression indicates that the 
loglognormal divergence can be observed in the vicinity of K m . Specifically, (jSJ) yields 
an asymptotic expression of F(k) near K m , 

exp{-A 2 ln 2 [-iA 1 ln[^/(2/ tm )]]} 



F K] 



A, 



5K\hx[5K/{2K m )]\ ' ^ 

where 5k = n m —K and < k//« m < 1. We find that our numerically obtained PDF best 
fits this asymptotic expression with the constants Ai = 1.7696, A 2 = 0.8342, A 3 = 1.0983. 
For more clear demonstration of the asymptotic behaviour of numerically obtained PDF, 
we employ a method which allows us to map the small neighborhood of K m into a 
wider range. Following the procedure applied to analyze the PDFs for CQNLSE and 
DCQNLSE M, we rewrite ® as 



In 



— SK\lll[SK/(2K m )}\F{K) 

A3 



A 2 ln 2 



•-Ailn[5/c/(2/c m )] 



We now define G{5k) and g(5n) as 

1 



G(5k) = - In 



A, 



and 



g(Sn) 



A 2 In 2 



6K\\n[6K/(2K m )]\F(K) 



--Ailn[5/t/(2/t r 



(7) 



(8) 



(9) 



Notice that if the numerically obtained PDF can be described by the rhs of (JBJ), the 
graph of G vs g is a straight line with a slope close to 1. By plugging the numerically 
obtained PDF data into F(k) in (jHJ), we calculate G(5k) and present the graph of G(5k) 
versus g{5n) in figure HI The graph demonstrates that our numerically obtained data 
lies on a straight line with a slope 1.04. This result concludes that the numerically 
obtained PDF of k exhibits a loglognormal divergence in the vicinity of K m . 



4. Conclusion 



We investigated numerically the evolution of flat-top solitary waves by the extended 
Korteweg-de Vries equation. Taking into account the disorder in linear growth/damping 
coefficient, one of the most common disorder forms in the nonlinear wave system, we 
showed that the PDF of k which characterizes the solitary wave amplitude, exhibits 
loglognormal divergence near the maximum value of k. We expect this loglognormal 




Figure 3. The probability density function of k at t = 10. 
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Figure 4. G(8k) vs g(<5/c) for the same parameters used in figure [3] The solid line is 
a linear fit with slope 1.04. The circles represent the numerical result. 



divergence in the vicinity of maximum possible value of k is mainly related to the fact 
that the solitary wave forms the table-top shape as the amplitude of wave approaches 
its maximum. The eKdV equation fundamentally differs from NLS type of equations 
considered earlier, however, it shows similar statistical behaviour for flat-top solitary 
waves. This phenomena suggests that the loglognormal divergence can be a general 
consequence associated with flat-top solitary waves in the presence of weak dissipative 
disorder. 
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